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We study the capture of a diffusing "lamb" by diffusing "lions" in one dimension. The capture 
dynamics is exactly soluble by probabilistic techniques when the number of lions is very small, and 
is tractable by extreme statistics considerations when the number of lions is very large. However, 
the exact solution for the general case of three or more lions is still not known. 



O I. INTRODUCTION 

CN . 

What is the survival probability of a diffusing lamb vifhich is hunted bw N hungry lions? Although this capture 
r— I process is appealingly simple to define (see Fig. |l|), its long-time behaviortTcl poses a theoretical challenge because 
1^ of the delicate interplay between the positions of the lamb and the closest lion. This model also illustrates a general 
^ feature of nonequilibrium statistical mechanics: life is richer in low dimensions. For spatial dimension d > 2, it is 
known that the capture process is "unsuccessful" (in the terminology of Ref. 1), as there is a nonzero probability for 
^ , the lamb to survive to infinite time for any initial spatial distribution of the lions. This result is a consequence of the 
transience of diffusion for d > 2,Dij which means that two nearby diffusing particles in an unbounded d > 2 domain 
may never meet. For d = 2, capture is "successful", as the lamb dies with certainty. However, diffusing lions in d = 2 
^ i are such poor predators that the average lifetime of the lamb is infinite! Also, the lions are essentially independent, El 
' so that the survival probability of a lamb in the presence of N lions in two dimensions is S'Ar(t) oc Si{t)^ , where Si{t), 
. the survival probability of a lamb in the presence of a single lion, decays aM (Int)^^. 
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' Lions are more efficient predators in d = 1 because of the recurrence of diffusion,HEI which means that two diffusing 

O • particles are certain to meet eventually. The d = 1 case is also special because there are two distinct generic cases. 
^ When the lamb is surrounded by lions, the survival probability at a fixed time decreases rapidly with N because 
the safe zone which remains unvisited by lions at fixed time shrinks rapidly in N . This article focuses on the more 
interesting situation of N lions all to one side of the lamb (Fig. |l|), for which the lamb survival probability decays as 
^ a power law in time with an exponent that grows only logarithmically in N . 

We begin by considering a lamb and a single stationary lion in Section IL The.|Survival probability of the lamb 
S\{i) is closely related to the first-passage probability of one-dimensional diffusionatl and leads to S\(t) ^ t~^^^. It 
is also instructive to consider general lion and lamb diffusivities. We treat this two-particle system by mapping it 
onto an effective single-particle diffusion problem in two dimensions with an absorbing boundary to account for the 
death of the lamb when it meets the lion, and then solving the two-dimensional problem by the image method. We 
apply this approach in Section III by mapping a diffusing lamb and two diffusing lions, onto a single diffusing particle 
within an absorbing wedge whose opening angle depends on the particle diffusivities, and then solving the diffusion 
problem in this absorbing wedge by classical methods. 



FIG. 1. Space-time evolution in one dimension of A'^ = 4 diffusing lions (dotted lines) which all start at a; = and a single 
diffusing lamb (dashed) which starts at x = xq. The trajectory of the closest ("last") lion, whose individual identity may 
change with time, is indicated by the heavy solid path. 



1 



In Section IV, we study N ^ 1 diffusing lionsfl'EI An essential feature of this system is that the motion of the 
closest ("last") lion to the lamb is biased towards the lamb, even though each lion diffuses isotropically. The many- 
particle system can. be recast as a two-particle system consisting of the lamb and an absorbing boundary which, from 
extreme statisticSjLJ moves to the right as \/ADi,t In TV, where D^is the lion diffusivity. Because this time deoeijjdence 
matches that of the lamb's diffusion, the survival probability depends intimately on these two motions ,Lrll2l with 
the result that SN{t) ^ t~^" and (3n oc IniV. The logarithmic dependence of (3^ on N reflects the fact that each 
additional lion poses a progressively smaller marginal peril to the lamb — it matters little whether the lamb is hunted 
by 99 or 100 lions. Amusingly, the value of (3n implies an infinite lamb lifetime for iV < 3 and a finite lifetime 
otherwise. In the terminology of Ref. 1, the capture process changes from "successful" to "complete" when iV > 4. 
We close with some suggestions for additional research on this topic. 

II. SURVIVAL IN THE PRESENCE OF ONE LION 

A. Stationary Lion and Diffusing Lamb 

We begin by treating a lamb which starts at xq > and a stationary lion at a; = 0. In the continuum limit, the 
probability density p{x, t) that the lamb is at any point a; > at time t satisfies the diffusion equation 

dp{x, t) d^pjx, t) 



where is the diffusivity (or diffusion coefficient). The probability density satisfies the boundary condition p{x — 
0, t) = to account for the death of the lamb if it reaches the lion at a; = 0_and the initial condition p{x,t = 0) = 
6{x — xq). Equation (|l|) may be easily solved by the familiar image method.B For x > 0, p{x,t) is the superposition 
of a Gaussian centered at xq and an "image" anti-Gaussian centered at — a;o- 



p{x,t) 



1 



-{x+xof /iDit 



(2) 



The image contribution ensures that the boundary condition at a; = is automatically satisfied, while the full solution 
satisfies both the initial condition and the diffusion equation. Thus Eq. (^ gives the probability density of the lamb 
for a; > in the presence of a stationary lion at a; = 0. 

The probability that the lamb meets the lion at time t equals the diffusive ffux to a; = at time t. The flux is 



dp{x,t) 
dx 



l/iDit 



(3) 



The ffux F(t) is also the first-passage probability to the origin, namely, the probability that a diffusing lamb which 
starts at xq reaches a; = for the first time at time t. Note that in the long time limit, defined by Dit ^ a;§, the 
first-passage probability reduces to F{t) x^jt^l"^ . This time dependence is a characteristic feature of the 

first-passage probability in one dimcnsion.u 

The probability that the lamb dies by time t is the time integral of F(t) up to time t. The survival probability 
is just the complementary fraction of these doomed lambs, that is. 



1- / F{t')dt', 
Jq 



^" (4) 
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The integral in Eq. (|J) can be reduced to a standard form by the substitution u = Xq/ y/ADit' to give 

Si{t) = erf f -^^] ~ as t oo, (5) 

where erf(z) — (2/^/%) e~" du is the error function.0 The same expression for Si{t) can be obtained by integrating 
the spatial probability distribution in Eq. (||) over all a; > 0. 

An amusing feature of the t~^^'^ decay of the lamb survival probability is that although the lamb dies with 
certainty, its average lifetime, defined as (t) = tF{t) dt — S{t)dt ~ t^^^'^ dt, is infinite. This infinite 
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lifetime arises because the small fraction of lambs which survive tend to move relatively far away from the lion. More 
precisely, the superposition of the Gaussian and anti-Gaussian in Eq. (|^) leads to a lamb probability distribution 
which is peaked at a distance {DitY/^ from the origin, while its spatial integral decays as {Dft)^^/"^ . 

B. Both Species Diffusing 

What is the survival probability of the lamb when the lion also diffuses? In the rest frame of the lamb, the 
lion now moves if either a lion or a lamb hopping event occurs, and their separation diffuses with diffusivity equal 
to + Dl (see for example, Ref. ||), where \s the lion diffusivity. From the discussion of Section IIA, the lamb 
survival probability has the asymptotic time dependence Si{t) ~ j ^ ixiD g + 0^)1. 
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FIG. 2. Mapping of the lion and lamb coordinates in one dimension to the planar coordinates j/i = xl/V Dl and 
2/2 = xe/\/De. The initial y-coordinates of the lion-lamb pair, (0, ^Di), and its image are indicated by the solid and open 
circles, respectively. Survival of the lamb, viVDl < y2\/De, translates to the diffusing particle in the plane remaining above 
and to the left of the absorbing line yw/Th] = y2\/Di. 



It is also instructive to determine the spatial probability distribution of the lamb. This distribution may be 
found conveniently by mapping the two-particle interacting system of lion at xl and lamb at xe, in one dimension 
to an effective single-particle system in two dimensions and then applying the image method to solve the latter (see 
Fig. H). To construct this mapping, we introduce the scaled coordinates yi = xl/-/IJl and t/2 = xi/y/Di to render the 
two-dimensional diffusive trajectory (2/1,2/2) isotropic. The probability density in the plane, p(yi,y2,t), must satisfy 
an absorbing boundary condition when y2VD£ — yi^/Dh, corresponding to the death of the lamb when it meets the 
lion. For simplicity and without loss of generality, we assume that the lion and lamb are initially at 2:^(0) = and 
xi{0) = 1 respectively, that is, 2/1 (0) = and 2/2(0) = \/TTi. The probability density is therefore the sum of a Gaussian 
centered at (2/1(0), 2/2(0)) = (0, ^/Dg) and an anti-Gaussian image. From the orientation of the absorbing boundary 
(Fig. 0), this image is centered at {\/DiSin26, —\/Di cos26'), where 9 ~ tanT^ D^/ D^. 




FIG. 3. Probability distribution of the lamb in one dimension at time t = IQ (Eq. (g)) when the lion and lamb are initially 
at a;_L = and xi = 1, respectively. The cases shown are r — Di/Dl = O.I, I, and 10 (bottom to top). 



From this image representation, the probability density in two dimensions is 



p{yi,V2,t) 



1 

47rt 



g-fe? + to2-\/57)^]/4t _ g-[(yi-\/D7sin2e)^+(j/2 + \/D7cos2e)2]/4t 



(6) 



The probability density for the lamb to be at 2/2 is the integral of the two-dimensional density over the accessible 
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range of the lion coordinate yi: 



p{y2,t) 



y2 cot 6 



(7) 



If we substitute the resuh (^) for p{yi,y2,t), the integral in Eq. c an be expressed in terms of the error function. 
We then transform back to the original lamb coordinate Xi — y2\/Di by using p{xe, t) dxt — p{y2, t) dy2 to obtain 



p{xe,t) 



1 



xi cot 9 \ 



) 



-(x£+cos 26)"^ /ADtt 



erfc 



'sin 29 — Xi cot 9 



(8) 



where erfc(z) = 1 — erf (z) is the complementary error function. A plot of p{xi, t) is shown in Fig. ^ for various values 
of the diffusivity ratio r = Di/Dl. The figure shows that the survival probability of the lamb rapidly decreases as 
the lion becomes more mobile. Note that when the lion is stationary, 9 = 0, and Eq. (^) reduces to Eq. (H). 

III. TWO LIONS 

To find the lamb survival probability in the presence of two diffusing lions, we generalize the above approach to 
map the three-particle interacting system in one dimension to an effective single diffusing particle in three dimensions 
with boundary conditions that reflect the death of the lamb whenever a lion is encountered.El Let us label the lions as 
particles 1 and 2, and the lamb as particle 3, with respective positions xi, X2, and X3, and respective diffusivities Di. 
It is again useful to introduce the scaled coordinates yi = Xij \fDi which renders the diffusion in the yi coordinates 
spatially isotropic. In terms of yi, lamb survival corresponds to 2/2V-D2 < ysV^s and yw/Di < ys^/D^. These 
constraints mean that the effective particle in three-space remains behind the plane 2/2\/^2 = 2/3V-D3 and to the 
left of the plane yiy/Di = ys^/D^ (Fig. ^(a)); this geometry is a wedge region of opening angle Q defined by the 
intersection of these two planes. If the particle hits one of the planes, then one of the lions has killed the lamb. 
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FIG. 4. (a) Mapping between the coordinates Xi of three diffusing particles on the line and a single isotropically diffusing 
particle in the three-dimensional space yi — Xi/\/TTi, subject to the constraints yi\/Di < ys^/Th, and y2\/D2 < ysVDa. The 
lamb survives if it remains within the wedge-shaped region of opening angle O. (b) Projection of the wedge onto a plane 
perpendicular to the e axis defined by the intersection of the two planes. 



This mapping therefore provides the lamb survival probability, sinccj-it is known that the survival probability of 
a diffusing particle within this absorbing wedge asymptotically decays asEJ 
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(9) 



For completeness, we derive this asymptotic behavior by mapping the diffusive system onto a corresponding electro- 
static system in Appendix A. To determine the value of Q which corresponds to our 3-particle system, notice that 
the unit normals to the planes yw/Ui = ys^/Ds and y2\/D2 — ysV^s are ni3 = (— 0, \/ \J D\ + D-^ and 
n23 = (0, — a/I?2, \fDl) I \/ D2 + fa, respectively. Consequently cos0 — fiia • (Fig. ^(b)), and the wedge angle is 
Q = 'K — (j) = -K — cos~^[D^/ [Di + 133) (-D2 + -D3)]. If we take Di = D2 = Dl for identical lions, and i^s = D^, the 
survival exponent for the lamb is 



/32(r) 



TT 

29 



2 cos"^ 



(10) 



where r = D^/Dl. 
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FIG. 5. The survival exponent (i2(r) given by Eq. (lid) versus the diffusivity ratio ' 



The dependence of P2{i') on the diffusivity ratio r is shown in Fig. ^. This exponent monotonically decreases from 
1 at r = to 1/2 for r — > 00. The former case corresponds to a stationary lamb, where the two lions are statistically 
independent and S2{t) = Si{t)'^. On the other hand, when r — + 00 the lamb diffuses rapidly and the motion of the 
lions becomes irrelevant. This limit therefore reduces to the diffusion of a lamb and a stationary absorber, for which 
S2{t) = Si{t). Finally, for De = Dl, P2 = 3/4 < 2/3i, and equivalently, S2{t) > Si{t)^. This inequality reflects the 
fact that the incremental threat to the lamb from the second lion is less than the first. 

IV. MANY LIONS 

The above construction can, in principle, be extended by recasting the survival of a lamb in the presence of TV lions 
as the survival of a diffusing particle in + 1 dimensions within an absorbing hyper- wedge defined by the intersection 
of the A^ half-spaces Xi < xn+i, * pf 1, 2, . . . , A. This approach has not led to a tractable analytical solution. On the 
other hand, numerical simulationsEl indicate that the exponent /3jv grows slowly with A", with (3^ « 0.91, P4 « 1.03, 
and PiQ w 1.4. The understanding of the slow dependence of Pn on A^ is the focus of this section. 

A. LOCATION OF THE LAST LION 

One way to understand the behavior of the survival probability is to focus on the lion closest to the lamb, because 
this last lion ultimately kills the lamb. As was shown in Fig. |l|, the individual identity of this last lion can change 
with time due to the crossing of different lion trajectories. In particular, crossings between the last lion and its left 
neighbor lead to a systematic rightward bias of the last lion. This bias is stronger for increasing A^, due to the larger 
number of crossings of the last lion, and this high crossing rate also leads to cciast (t) becoming smoother as A^ increases 
(Fig. |6|) . This approach of the last lion to the lamb is the mechanism which leads to the survival probability of the 
lamb decaying as t~^" , with f3]\[ a slowly increasing function of A^. 
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FIG. 6. Time dependence of xiast for a single realization of A'' = 4, 64, and 1028 lions (bottom to top). This data was 
generated by tracking the position of the rightmost among A'^ lions, each of which performs a nearest-neighbor discrete-time 
random walk starting from x — 0. The coincidence of the data and the linear early-time growth of xist are artifacts of the 
discrete random walk motion. The straight line of slope 1/2 indicates the expected long-time behavior. 



To determine the properties of this last lion, suppose that ^ 1 lions are initially at the origin. If the lions 
perform nearest-neighbor, discrete-time random walks, then at short times, xiast(i) = t. This trivial dependence 
persists as long as the number of lions at the last site in their spatial distribution is much greater than one. In this 
case there is a large probability that one of these lions will hop to the right, thus maintaining the deterministic growth 
of xiast- This growth will continue as long as --^===e~* /^^i-* ^ l that is, for t <C 4£'iln7V. At long times, an 
estimate for the location of the last lion is provided by the condition!] 



Equation ( pi] ) specifics that there is one lion out of an initial group of lions which is in the range [xiast , oo] . Although 
the integral in Eq. (^) can be expressed in terms of the complementary error function, it is instructive to evaluate it 
explicitly by writing x = xiast + e and re-expressing the integrand in terms of e. We find that 



Over the range of e for which the second exponential factor is non-negligible, the third exponential factor is nearly 
equal to unity. The integral in Eq. ([l^ ) thus reduces to an elementary form, with the result 

^^e---/4^-*^ = l. (13) 



If we define y = xust/V^^^U and M — N/^/Att, the condition in Eq. ( p^ can be simplified to 

ye^' = M, (14) 



with the solution 



lln(lnM) 



y = VlnM { 1 

^ V 4 InAf 



(15) 



In addition to obtaining the mean location of the last lion, extreme statistics can be used to find the spatial probability 
of the last lion. For completeness, this calculation is presented in Appendix B. 
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To lowest order, Eq. (|T^) gives 



xi^stit) ^ V4jJiMNt = yi^, (16) 



for finite N. For N = oo, xiast(0 would always equal t if an infinite number of discrete random walk lions were 
initially at the origin. A more suitable initial condition therefore is a concentration cq of lions uniformly distributed 
from — oo to 0. In this case, only N oc y^c^DiJ of the lions are "dangerous," that is, within a diffusion distance 
from the edge of the pack and thus potential candidates for killing the lamb. Consequently, for N ^ oo, the leading 
behavior of a;iast (t) becomes 



a:iast(t) - ^2Dl H4DLt) t . (17) 

As we discuss in Section IVB, the survival probability of the lamb in the presence of many lions is essentially determined 
by this behavior of xiast- 

B. LAMB SURVIVAL PROBABILITY FOR LARGE N 

An important feature of the time dependence of xiast is that fluctuations decrease for large N (Fig. |^). Therefore 
the lamb and N diffusing lions can be recast as a two-body system of a lamb and an absorbing boundary which 
deterministically advances toward the lamb as xust(t) = ^/AnI. 

To solve this two-body problem, it is convenient to change coordinates from [x, t] to [x' — x — xig^tit), t] to fix the 
absorbing boundary at the origin. By this construction, the diffusion equation for the lamb probability distribution 
is transformed to the convection-diffusion equation 

dp{x',t) _x^ dp{x' t) ^ M^' t) ^ (0<x'<oo) (18) 
ot 2t ox ax^ 

with the absorbing boundary condition p(x' = 0, t) = 0. In this reference frame which is fixed on the average position 
of the last lion, the second term in Eq. (|8|) accounts for the bias of the lamb towards the absorber with a "velocity" 
— a;iast/2t. Because xiast ^ V^Nt and x' ~ ^/Dit) have the same time|-jdeBendence, the lamb survival probability 
acquires a nontrivial dependence on the dimensionless parameter ^jv/ZJ^H IIJ Such a dependence is in contrast to the 
cases xiast ^ x' or xiast 3> x\ where the asymptotic time dependence of the lamb survival is controlled by the faster 
of these two co-ordinates. Such a phenomenon occurs whenever there is a coincidence of fundamental length scales in 
the system (see, for example, Ref. |T3[ ). 

Equation (|l^) can be transformed into the parabolic cylinder equation by the following steps.i First introduce 
the dimensionless length ^ = x'/xiast and make the following scaling ansatz for the lamb probability density, 

p(a;',i) -t-^"-i/2^(0. (19) 

The power law prefactor in Eq. (^) ensures that the integral of p(x', t) over all space, namely the survival probability, 
decays as t~^'^ , and F{^) expresses the spatial dependence of the lamb probability distribution in scaled length units. 
This ansatz codifies the fact that the probability density is not a function of x' and t separately, but is a function 
only of the dimensionless ratio x'/xiast{t)- The scaling ansatz provides a simple but powerful approach for reducing 
the complexity of a wide class of systems with a divergent characteristic length scale as t ^ ooJl3 



If we substitute Eq. (|19| ) into Eq. (18), we obtain 



Now introduce ri = + 1)^A^/2D^ and F{^) = e"'''/'* ^(t;) in Eq. This substitution leads to the parabolic 

cylinder equation of order 2/3ArliJ 



dry 2 



1 ry 



2?2/3„ = 0, (21) 



subject to the boundary condition, 2?2/3]v('?) — fo^' both ?/ = y/ A^/^Di and ry = oo. Equation (|2^) has the form 
of a Schrodinger equation for a quantum particle of energy 2/3 jv + ^ in a harmonic oscillator potential 77^/4 for 
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rj > ^ Ajq l2Di^ but with an infinite barrier at 77 = ^ Aj<i jDix^ For the long time behavior, we want the ground state 
energy in this potential. For ^ 1, we may approximate this energy as the potential at the classical turning point, 
that is, + 5 — '7^/4. We therefore obtain /3jv ^ jX^Dg. Using the value of Aj^ given in Eqs. ( |l6| ) and ( [T7| ) 
gives the decay exponent 

f In TV, if finite 
Hn^K ' (22) 
[^ln<. foriV = oo 



The latter dependence of (i-^ implies that for N — > 00, the survival probability has the log-normal form 



5ooW^exp( In^i) . (23) 



Although we obtained the survival probability exponent /3jv for arbitrary diffusivity ratio r = Di/Dl, simple 
considerations give different behavior for r 3> 1 and r <C 1. For example, for r — Q (stationary lamb) the survival 
probability decays as t~^/'^. Therefore Eq. ( p^ ) can no longer apply for r < N^'^, where /3Ar(r) becomes of order 
N. Conversely, for r — 00 (stationary lions), the survival probability of the lamb decays as t~^/^. Thus Eq— ([2^) 
will again cease to be valid for r > In A^, where /3Ar(r) becomes of order unity. By accounting for these limits,B the 
dependence of on the diffusivity ratio r is expected to be 

f N/2 r<f/7V 
/3^(r) = <^ (l/4r)ln(4iV) 1/N r \n N (24) 

The r dependence of (3n in the intermediate regime of 1/A^ ^ r <C InA^ generalizes the exponents given in Eq. ( |l0| ) 
for the three-particle system to general N. 

V. DISCUSSION 

We investigated diffusive capture of a lamb in one dimension by using of several essential techniques of nonequi- 
librium statistical mechanics including first-passage properties of diffusion, extreme value statistics, electrostatic 
analogies, scaling analysis, and moving boundary value theory. These tools provide an appealing physical description 
for the survival probability of the lamb in the presence of N lions for the cases of very small and very large N. 
Nevertheless, the exact solution to capture of the lamb remains elusive for iV > 3. 

We close by suggesting several avenues for further study: 

1 . Better Simulation Methods. Previous simulations of this systemB followed the random walk motion of one lamb 
and N lions until the lamb was killed. This type of simulation is simple to construct. One merely places 
the lamb and lions on a one-dimensional lattice and have them perform independent nearest-neighbor random 
walks until one lion lands on the same site as the lamb. The survival probability is obtained by averaging 
over a large number of realizations of this process. However, following the motion of discrete random walks is 
inefficient, because it is unlikely for the lamb to survive to long times and many realizations of the process must 
be simulated to obtain accurate long-time data. In principle, a_much better approach would be to propagate 
the exact probability distribution of the particles in the systemJlE Can such an approach be developed for the 
lamb-lion capture process? Another possibility is to devise a simple discrete random-walk process to simulate 
the motion of the last lion. Such a construction would permit consideration of just the lamb and the last lion, 
thus providing significant computational efficiency. 

2. The Last Lion. Extreme value statistics provides the spatial probability distribution of the last lion. We may 
also ask other basic questions: How long is a given lion the "last" one? How many lead changes of the last lion 
occur up to time tl How many different lions may be in the lead up to time t? What is the probability that a 
particular lion is never in-the lead? Methods to investigate some of these issues are also outlined in the article 
by Schmittmann and ZiallZl in this journal issue. 

3. Spatial Probability Distribution of the Lamb. As we have seen for the case of one lion, the spatial distribution 
of the lamb is a useful characteristic of the capture process. What happens for large iV? In principle, this 
information is contained in the solution to the parabolic cylinder equation for the scaled probability distribution 
(Eq. (|2^)). The most interesting behavior is the form of the distribution close to the absorbing boundary, where 



8 



the interaction between the lions and the lamb is strongest. For N = 1 lion, this distribution decays linearly to 
zero as a function of the distance to the lion, while for N = 2, the distribution has a power law decay in the 
distance to the last lion which depends on the diffusivity ratio Di/Dl. What happens for general N and for 
general diffusivity ratio? Is there a physical way to determine this behavior? 

4. Two-Sided Problem. If N lions are located on both sides of the lamb, then the lamb is relatively short-lived 
because there is no escape route. One can again construct a mapping between the N + 1-particle reacting 
system and the diffusion of an effective particle in an absorbing wedge-shaped domain in + 1 dimensions. 
From this mapping, S'Ar(t) decays as t~'^^ , but the dependence of 7Ar on TV and diffusivity ratio is unknown. It 
is clear, however, that the optimal strategy for the surrounded lamb is to remain still, in which case the lions 
are statistically independent and we then recover SN{t) ^ t^^^"^ . Is there a simple approach that provides the 
dependence of 7j\f 'ijK.A^ for arbitrary diffusivity ratio? Finally, Sao{t) exhibits a stretched exponential decay in 
time, exp(— What is that nature of the transition from finite N to infinite N behavior? 

5. Intelligent Predators and Prey. In a more realistic capture process, lions would chase the lamb, while the lamb 
would attempt to run away. What are physically-reasonable and analytically-tractable rules for such directed 
motion which would lead to new and interesting kinetic behaviors? 
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APPENDIX A: SURVIVAL PROBABILITY OF A DIFFUSING PARTICLE WITHIN A WEDGE 

The survival probabilitK-pf a diffusing particle within an absorbing wedge can be derived by solving the diffusion 
equation in this geometryB-EJ We provide an alternative derivation of this result can also be obtained by developing 
a correspondence between diffusion and electrostatics in the same boundary geometry. Although the logic underlying 
the correspondence is subtle, the result is simple and has wide applicability. 

The correspondence rests on the fact that the integral of the diffusion equation over all time reduces to the 
Poisson equation. This time integral is 

m-v[rAt)-'-^^]dt. (25) 



If one defines an electrostatic potential by (f>(r, 8) — p{r, 9, t) dt, Eq. (25) can be written as 



VHir, 0)^-^ b(r, 9,t = ^)- p{r, 0, t ^ 0)] . (26) 

For a boundary geometry such that the asymptotic survival probability in the diffusive system is zero, then Eq. ( p6| ) is 
just the Poisson equation, with the initial condition in the diffusive system corresponding to the charge distribution in 
the electrostatic system, and with absorbing boundaries in the diffusive system corresponding to grounded conductors. 

To exploit this analogy, we first note that the electrostatic potential in the wedge decays as <I>(r, 6) ^ r^^/^ for 
r — > oo, for any localized charge distribution.E3 Because the survival probability of a diffusing particle in the absorbing 
wedge is given by S{t) = J p{r, 6, t) dA, where the integral is over the area of the wedge, we find the following basic 
relation between S{t) and the electrostatic potential in the same boundary geometry 

S{t)dt ^ [ dt [ p{r,0,t)dA 



"'0 



/Dt /.e 







rdr dO $(r, 
Jo 



/Dt 







oct^-''/^^. (27) 

In evaluating the time integral of the survival probability, we use the fact that particles have time to diffuse to radial 
distance ^/Dt but no further. Thus in the second line of Eq. (p7|), the time integral of the probability distribution 
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reduces to the electrostatic potential for r < \/ Dt and is essentially zero for r > V Dt. Finally, by differentiating the 
last equality in Eq. (|2^) with respect to time we recover Eq. (^). 

APPENDIX B: SPATIAL DISTRIBUTION OF THE LAST LION BY EXTREME STATISTICS 

It is instructive to apply extreme statistics to determine the probability distribution for the location of the last 



lion from an ensemble of N lionsE Let p{x) = 



g X /ADLt -j-j^g (^Gaussian) probability distribution of a 



single lion. Then py{x) = p(x')dx' is the probability that a diffusing lion is in the range [a;,(X)] and similarly 
p<{x) = 1 — p-^(x) is the probability that the lion is in the range [—00,0:]. Let L]sf{x) be the probability that the last 
lion out of a group of TV is located at x. This extremal probability is given by 



Ln{x) = Np{x)p<{x) 



(28) 



That is, one of the N lions is at a;, while the remaining 1 lions are in the range [— 00, x]. If we evaluate the factors 
in Eq. (Eq), we obtain a double exponential distribution:EI 



Ln(x) 



N 



N 



-x^/ADlI 



-x-'/4DLt 



exp 



1 

N ~l 



N-l 



^/A^njrtjx 



(29) 



When TV is large, then x/ \/ADi,t is also large, and we can asymptotically evaluate the integral in the exponential in 
Eq. (^9|). Following Eq. (|lj), it is convenient to express the probability distribution in terms of M — Nj^f^ and 
V — a;iast /\/4i7tI. If we use Lf^(%)) dy — Ljs[{x) dx, we obtain 



ijv(y) ~ 2Me-^ exp(-Me-^ /y). 



(30) 



The most probable value of xiast is determined by the requirement that L'j^{y) = 0. This condition reproduces 



M given in Eq. (|14|) 



We may also estimate the width of the distribution from its inflection points, that is, when L'^j{y) = 0. By 



straightforward calculation, L'^{y) = at 



y± 



y/\n{M/k±) w VlnM(^l 



In k± 
InM 



(31) 



where k± — (3± V5)/2. Therefore as ^ 00, the width of LN{y) vanishes as l/VhiM. This behavior is qualitatively 
illustrated in Fig. ^ where the fluctuations in xiast(0 decrease dramatically as N increases. This decrease can also be 
understood from the form of the extreme distribution Ln{x) in Eq. (|29|). The large-x decay of Ln{x) is governed by 
p{x), while the double exponential factor becomes an increasingly sharp step at a;stop ~ ^/^D^t In as N increases. 
The product of these two factors leads to L^^x) essentially coinciding with p{x) for x > Xstop and Lj^{x) « for 



X < X. 
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